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EXCITED-STATE ENERGIES AND SUPERSYMMETRIC 

INDICES 

PAUL FENDLEY 


Abstract. We study multi-soliton states in two-dimensional N =2 supersym¬ 
metric theories. We calculate their energy exactly as a function of mass 
and volume in the simplest integrable N=2 supersymmetric theory, the sine- 
Gordon model at a particular coupling. These energies are related to the 
expectation value X = tr[exp(m7ri 7 ') exp (-H/T)], where F is the fermion num¬ 
ber. For n=l, this is Witten’s index; for n an odd integer, we argue that X is 
an index in the sense that it is independent of all D-term variations. 


1. Introduction 


The thermodynamic Bethe ansatz (TBA) has proven a useful tool for computing 
the ground-state (Casimir) energy of an integrable 1 + 1 dimensional system |Q. 
The special properties implied by integrability allow one to compute the partition 
function at any temperature T 


(i) 


Z = tr 


e ~H L /T 


where Hr is the Hamiltonian for the system on a line of length L. The TBA 
calculation gives the free energy per unit length —T\nZ/L when L much is larger 
than any other length scale in the theory. This partition function is related to the 
ground-state energy by interchanging the definition of space and time. The above 
trace is equivalent to a Euclidean path integral where space time is a cylinder of 
length L and circumference 1/T. This path integral can equivalently by computed 
by taking space to be a circle of circumference R = 1/T, and time propagation to 
be in the L-direction: 

(2) Z = tr [e~ H * L ] 


Since L is very large, the leading contribution to this trace is from the lowest 
eigenvalue of Hr, which is the ground-state energy Eq. Thus 

(3) E 0 = -^\nZ 

This trick of reversing space and time is known in the conformal field theory liter¬ 
ature as a modular transformation. 

The TBA can be generalized in order to calculate excited-state energies as well. 
Adding imaginary chemical potentials to ( 0 ) corresponds to changing the boundary 
conditions around the cylinder (i.e. putting some sort of twist in the R-direction). 
This enables one to project out the ground state [f| and calculate excited-state 
energies [||]. In this manner one can calculate a few low-lying excited-state energies 
in a number of models [Q, [|, but all those calculated here have the property that 
for large R, they are degenerate with the ground state and thus not multi-particle 
states. 
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Recently, several techniques were developed which allow the calculation of the 
energies of multi-particle states §§■ The basic idea is to analytically continue the 
TBA equations onto another sheet in an appropriate space (the fugacity in Q , the 
mass in ||). Using formal techniques developed in Q, the results of the continuation 
can be written in closed from, and are interpreted as the TBA equations for an 
excited state. These generalized TBA equations have been worked out in detail 
only for several states in the Yang-Lee model || |[ and some generalizations ||. 

We give explicit thermodynamic Bethe ansatz equations for excited-state energies 
in a unitary theory, the simplest integrable N=2 supersymmetric field theory. In 
this model, which is equivalent to the ordinary sine-Gordon model at a special value 
of the coupling, we calculate an infinite number of eigenvalues of the Hamiltonian, 
namely 

(4) E n = (n\H R \n) 

when n is a multi-soliton states on a circle of circumference R. The energy E n 
of this excited state is a function of the volume R of the one-dimensional space 
and the mass scale m of the theory. The ground-state energy Eq is the Casimir 
energy on this circle, and it vanishes as R —> oo. The excited states in general are 
a collection of interacting particles, so as P —> oo, E n —+ am, with a depending on 
how many particles are in the state and their masses. If there is only one particle 
mass in the theory, then a is an integer. 

At first glance, the supersymmetric model to be discussed here is more compli¬ 
cated than the Yang-Lee model studied in detail in §!• It has two particles in the 
spectrum (the soliton and antisoliton in sine-Gordon) with a non-diagonal S'-matrix 
instead of a single particle with a single S'-matrix element. However, because of 
various simplifying features probably related to the supersymmetry, we are able to 
find the appropriate equations to describe an infinite number of excited-state ener¬ 
gies. The excitations in this theory are all solitons saturating the Bogomolny bound 
Jq| , so our results give the exact non-perturbative energy of an infinite hierarchy of 
interacting multi-soliton states. 

Some of these excited-state energies have a particular interest in the context of 
N= 2 theories. Formally, one can view the calculation in terms of 

(5) l(a) = ti[e iaF e~ RHL ], 

corresponding to adding imaginary chemical potentials proportional to the fermion 
number F. T(a) makes sense as a function of a at least for a small, and is real 
when the theory is CP-invariant. It depends on R, L and any mass scale m in the 
theory. By analogy with (||) we define 

(6) lnX(a) = —E(a)L 

for large L. Below we discuss how to define E(a) for all a. The point of |fj is that 
at appropriately-chosen values a n , one has 

(7) E n = E(a n ) 

In this supersymmetric case, we study a n = mr. We show that when appropriately 
defined, 1(a) is not periodic in a , even though F is half-integer in our example. 
(This arises because of subtleties in taking the infinite-space limit, as are common 
in open-space index theorems JlTij] .) 
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In supersymmetric theories, an object of fundamental importance is Witten’s 
index @ 


( 8 ) 


1(tt) = tr 


(_1 ) F e ~ HL /T 


In the case discussed here, Witten’s index is independent of R , L and m. Another 
case of interest in N=2 theories is d/(a)/da| a=7r , which is the “new supersymmetric 
index” calculated in |l2) . This is an index in the sense that it is independent of all 
-D-terrn perturbations of the theory. In this paper, we will show that I(mr) is also 
independent of the H-term when n is an an odd integer. 

More precisely, in section 2 we show that I(nir) for n an odd integer corre¬ 
sponds to computing Witten’s index with different open-space asymptotic condi¬ 
tions. These asymptotic conditions correspond to inserting operators at the ends 
of the space-time cylinder. We argue that these objects are independent of D-term 
type perturbations, and thus are open-space indices in the sense of Q. We also 
show the 1(a) gives excited-state energies when n is any integer. In section 3 we 
use the thermodynamic Bethe ansatz to derive equations for 1(a) in our example. 
In section 4, we discuss the solutions in the m —> 0 limit, where 1(a) is an analytic 
function of a. In section 5 we extend the results to all mR, showing how to define 
T(a) past discontinuities at certain values of a. In the first appendix, we check our 
results by explicitly computing the energy in the massless limit. In the second ap¬ 
pendix, we calculate the excited-state energies of a Dirac fermion, mainly in order 
to clarify the nature of the discontinuity arising in section 5. 


2. Supersymmetric indices and excited states 

In this section we explain how the excited-state energies discussed in the intro¬ 
duction are related to interesting quantities in N—2 supersymmetric field theories. 

In these two-dimensional models, Euclidean spacetime is a cylinder of circumfer¬ 
ence R and very long length L. Varying a in (|s]) corresponds to changing boundary 
conditions as one goes around the cylinder. In particular, for a = mr, when n is 
an even integer these correspond to antiperiodic boundary conditions on fermionic 
operators. When n is an odd integer, these are periodic boundary conditions on the 
fermions. In conformal held theory (which applies when m —> 0), these are called 
Neveu-Schwarz and Ramond boundary conditions respectively. 

The energy E(a) in (Q) is a natural object to consider when one thinks of space 
as in the indirection, and time in the L-direction. It is a particular eigenvalue of the 
Hamiltonian Hr operating on the space of states | ) a on a circle of circumference R. 
The Hilbert space is labeled by a because changing the boundary conditions around 
R changes the states allowed in the Hilbert space. Shifting a by 27 t takes the space 
of states | ) a to itself, so for example the Ramond sector is mapped to the Ramond 
sector. However, it does not take a given state to itself. This is well known in the 
superconformal literature where the change in a is often called spectral flow; see 
e.g. Q. Thus in the m —> 0 limit, E n for n odd is the energy of some state in the 
Ramond sector, while for n even it is the energy of some state in the Neveu-Schwarz 
sector. 

In a conformally-invariant theory, the energy of a state must be proportional to 
1 /R because there is no other scale in the theory. The constant of proportionality 
is related to the central charge c of the theory and the left and right conformal 
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dimensions (h n , h n ) of the operator which creates the state | n) Q: 

(9) \imE n = --^-(c-12h n -12h n ). 

m —>0 bK 

Thus in the m —* 0 limit, the ground-state energy in the conformal limit is related 
to the central charge, while excited-state energies are related to the dimensions of 
the operators which create these states. Thus when n is an even (odd) integer, 
we can identify which operator in the Neveu-Schwarz (Ramond) sector created the 
state. 

In an N=2 theory, spectral flow arguments give |li| 

(10) =-Si 0 - © 1 2 ) ■ 

Since calculations here are left-right symmetric, the states | n) are therefore created 
by operators with conformal dimension (h n ,h n ) = (cn 2 /24,cn 2 /24). The formula 
© is well known in the example we will discuss in later sections, the sine-Gordon 
model at its N=2 supersymmetric point. The m —> 0 limit consists of a free boson, 
which has c = 1. [] 

Notice that at n = 1, we have E\ = 0. In fact, E\ = 0 for any m, not just in 
the conformal limit (see [|l2| for the TBA calculation in many N=2 theories). This 
is a deep result for supersymmetric theories, because from the point of view where 
space is in the L-direction and time is in the .R-direction, we are calculating the 
leading term in Witten’s index [jTlJ] . Witten’s result is that if there are periodic 
boundary conditions in the L-direction, then X(j r) in (|s|) is related to the number 
of ground states, and is an integer independent of L,T and to. Thus when we find 
E by interchanging the roles of space and time, ®) shows that E\ = 0 because X is 
finite and L —* oo. With periodic boundary conditions in the L-direction, Witten’s 
index is independent of any supersymmetry-preserving deformations of the theory. 

In general X(a) does depend on perturbations, for example, changing the mass. 
However, when n is odd, X(mr) does have a special property. We argue here 
that it is independent of any D -term perturbations. There are four supercharges 
Q + ,Q~,Q + ,Q~ in an N=2 theory (for any to). Both Q + and Q have fermion 
number +1, while the other two have fermion number —1. These charges are de¬ 
fined with periodic boundary conditions, so for the rest of this section we restrict to 
the case n an odd integer. In a superspace action a D-terrn comes from integrations 
over all four Grassman coordinates. This means that a variation of a D-term can 
be written as inserting {Q + , [Q , A(x)]} in the path integral where A itself can be 
written as {Q ~, [Q , A']} for some K. Thus a H-ternr variation of X(mr) can be 
written as 

(11) SX^nn) = — Jdx(n\{Q + ,[Q , A(x)]}e~ HRL \n) 

Because Q + and Q anticommute with each other, if either Q + \n) = 0 or Q \n) = 
0, then SI = 0. 

In a superconformal field theory, a state annihilated by Q+ is taken to another 
one annihilated by Q + under spectral flow by positive a. We prove this by writ¬ 
ing the left-moving and right-moving fermion currents J(x + t) and J(x — t) in 

1 Conformal field theory aficionados will notice that (|uj) looks like the central charge of the 

minimal models. A non-zero a corresponds to a non-zero charge at infinity in the Coulomb-gas 
picture & Taking a = n gives the trivial c = 0 minimal model. 
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terms of a free boson: J = S x <$>(x,t), and J = — <5 x <I>(:r, t). The fermion number 
F = J^'dx[J(x +1) — J(x — f)]. In the sine-Gordon example to be discussed be¬ 
low, $ — $ is just the c = 1 boson, normalized so that the perturbing operator of 
dimension (2/3, 2/3) is cos[2(<h — $)]. Then Q + = J dxAexp(i3&(x,t)/c), while 
Q = J dxB exp(i3$(x, t)/c) 1 where A and B commute with $ and *I>. The ex¬ 
ponents can be verified by checking that they give the supercharges the correct 
fermion number. The operator which implements the spectral flow can be written 
in terms of this boson. The state | n) in the Ramond sector is [[Hi] 

r R 

| n)= / dx e is G > ( x > t )- 4> ( a: > i ))|i)_ 

Jo 

where s = {n — l)/2 is the number of “units” of spectral flow; one unit of spectral 
flow corresponds to shifting a by 27 t. Therefore, if s is positive, the spectral flow 
operator commutes with Q + , because the exponents of both are positive. If s is 
negative, the spectral flow commutes with Q . Since the state | n = 1) yields 
Witten’s index, it is annihilated by Q + and Q (and the other supercharges as 
well). Therefore, the states with n an odd integer are all annihilated by either Q + 
or Q . This argument applies only at the superconformal point. To go off the 
critical point while preserving supersymmetry, one adds to the Hamiltonian a term 
which commutes with all the supercharges. Therefore, in particular, it commutes 
with Q + , and 51 in ([ll]) is still zero. 

We have therefore shown that the quantities I(mr) are independent of any ZD- 
term variation of the theories when n is an odd integer. Generically, there are 
only a finite number of perturbations of a theory which are not ZD-term, so this 
means I(mr) depends on only a finite number of parameters. In the example we 
will discuss below, there is only one. 

One can think of l(nir) with n odd as an open-space index. Before Witten’s 
work, in fact, it was shown that there are theories where (^j) is not an integer and 
does depend on various parameters The arguments of 0 utilize periodic 

conditions and hence do not necessarily apply in open space. The expectation 
value (|^) can vary when defined on an open space with appropriate asymptotic 
conditions. We have shown that l(mr) with n odd amounts to calculating (||) 
with asymptotic conditions corresponding to inserting the state |n) at the end of 
the cylinder. For the subsequent analysis, we do not need to know how to write 
these boundary conditions in terms of the fields in Hl. Knowing how would be 
interesting, in order to compare our results with Jlo| and subsequent open-space 
index results like 0- 

Z(a) as defined in (|^) seems periodic in a when F is rational. Periodicity in a 
contradicts for example the result (0). Thus the definition (||) should be under¬ 
stood as the appropriate continuation away from a = 0. In particular, below we 
will calculate E(a ) for all a and use this to define X(a). As with the open-space 
index theorems, this seems to be the only way to precisely define I{a) for all a 
when L —> 00 . 


3. The thermodynamic Bethe ansatz equations 

In the previous section we showed that E(a) = — In 1(a)/L gives excited-state 
energies in N=2 supersymmetric theories when a = nir. For the remainder of this 
paper, we show how to calculate E n = E(nn) for the sine-Gordon model at its 
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N = 2 supersymmetric point. This is the first model in the N=2 discrete series 
perturbed by its only relevant supersymmetry-preserving operator (see e.g. [9) for 
a discussion of these models). This model is interesting not just because of its 
supersymmetry. Because it corresponds to sine-Gordon coupling /3 2 = 167r/3 (in 
the usual normalization where the kinetic term in the action has a 1/2 in front, 
and the free fermion point coupling /3 2 = 47 t), it is in the “repulsive” regime, where 
the fermions in the related Thirring model repel each other. The cases discussed in 
|5| |(i|, U are all in the attractive regime. Moreover, as we mention in the conclusion, 
our model gives information on 2d circular polymers. 

Our starting point is the thermodynamic Bethe ansatz (TBA) equations in func¬ 
tional form. Here, these equations are written in terms of two entire functions 
X(9; a) and Y(9; a) 

(12) X(9 + m/2)X(9 -iit/2) = l + 2cos aY(9) + Y 2 {9) 

(13) Y{0 + iir/2)Y(0 -in/2) = 1+X(0) 

where for ease of notation we often don’t write out the a argument. For a < 7r 
these equations can be derived from the underlying lattice model mu or directly 
from the quasiparticle S matrix || |l2|. It follows from ([l|) that when X ^ 0, 
Y(0) = Y{9 + 3*7t) and X{9) = X(9 + 3in). Even though the equations are periodic 
in a, there is no guarantee that the solutions will be. 

To relate these functions to the energy levels of the system, we utilize a conjec¬ 
ture of ||. This conjecture relates the TBA quantities to the infinite number of 
conserved quantities I n (a) in an integrable theory. In this case (/3 2 = 2/3 in their 
conventions), their conjecture relates the asymptotic expansion of Y(0\a) around 
9 = —oo to the I n : 

OO 

(14) Iny(0; a) = ]T C n e {2n ~ 1)e I 2n -i(a) 

n =1 

where the the convention-dependent constants C n are proportional to i? 2n_1 . The 
first conserved quantity Ii{a) is just the energy E(a). As we will see momentarily, 
this is equal to the ground-state energy determined by standard TBA methods. 

To make further contact with the physics, we note that 9 is the rapidity of the 
quasiparticles of the problem, defined so that when the quasiparticles have mass 
m, the energy of the particle is m cosh 9 and the momentum is to sinh 9. The filling 
fraction (the Fermi distribution function generalized to this interacting problem) at 
temperature T = 1/R is X/(l+X) when a < 7r [[)]. This means that the probability 
that there is a particle with rapidity between 9 and 9 + d9 is d9X(9)/( 1 + X(9)). 
At high energy (large 9) the particle interactions are irrelevant, so we can treat 
these particles as free and having a Fermi distribution function. This implies an 
asymptotic condition for X (9 ): 

(15) X(9) -f e ~ mRcoshe f or Q ±oo 
which in turn requires 

(16) Y(9) —> 1 for Re0 —> ±oo, \Im9\ < n/2 

Even with the asymptotic condition (|l5|), there are an infinite number of solutions 
to ©• To fix a unique solution, we utilize a lemma of |Q : 


EXCITED-STATE ENERGIES AND SUPERSYMMETRIC INDICES 


7 


Lemma: When f(9) is analytic and bounded in the strip —7 t/2 < Im(9) < 7t/ 2 
and obeys the relation 


f(6 + m/2) + f{9 - iir/2) = g(9) 


for some g(9), then 


(17) 


m = 



dff 1 
27t cosh (6 — 9') 


9 ( 0 ') 


Another derivation of this lemma is given in pdf . 

If there are no zeroes in X(0) or Y(0) in \Im0\ < n/2 and X obeys the asymptotic 
condition dH), then both ln[X(0)e m ' Rcoshe ] and In Y(d) satisfy the conditions of the 
lemma above. Thus there is a unique solution Xq(9),Yo(9) of © with no zeros 
and which satisfies the asymptotic condition. The lemma requires that they solve 
the integral equations: 


In A 0 (9) 

(18) 

(19) lnF o (0) 


—mR cosh 9 

+/ € within - I (1 ++ 

/ 2» cosh(0 - »') hi:1 + ' '' 


where all rapidity integrals here and for the rest of the paper run from — oo to 
oo unless otherwise labeled. These are the usual TBA equations for this system 
|l8l [l9| . P, so we conclude that for |a| < 7r there indeed are no zeroes in X(9) or 
Y(9) for \Im(9)\ < n/2. From © and (^9[) it follows that the ground-state energy 
is 


( 20 ) 


Ec) = — m 


d9 

2ir 


coshd ln(l + Xq(9)), 


the usual TBA form. 

As we will discuss in detail in the next section, X and Y develop zeroes when 
a > 7T. Thus the excited-state energies are given in terms of X and Y obeying more 
complicated TBA equations. To find these, we first note that 



Note also that in the region \Im9\ < 7t/ 2, tanh((0 — x)/2) is bounded and its only 
zero occurs at 9 = x. Thus when X(9) has zeroes in this region at 9 = xi, X 2 ,..., xj 
(and at 9 = —xi, — X 2 ,..., — xj due to the symmetry under 9 —> —9) and Y(9) has 
zeroes at 9 = y i,y 2 ,..., Yk and their opposites, then 


and 


In 


X{9)e mRcoshe coth 

3 =1 


9 — x, \ / 9 + x, 

1 coth 


In 


K 


y(9) n coth 


k =1 


0 -yk 


coth 


(^) 
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satisfy the conditions of the Lemma. Thus the solution of (|13|) for a given J and 
I\ and asymptotic form (lla) obeys the integral equations 


( 21 ) 


lnX($) = —mR cosh 0 + In 

j =i 

dO' 1 


tanh 


tanh 


( 


V 2 


+ 


27t cosh(0 — O') 


In [(1 + e ia Y(0’))( 1 + e~ ia Y(0'))] 


K 


lnF(6») = ln 


k—1 


( 22 ) 


tanh 

dff 


0-J k 
2 

1 


tanh 


0 + y k 


ln(l + X(0')) 


J 27t cosli(0 — O') 

The energy associated with this solution is then 

1 f dO ^ 

(23) E(a) = — — lnX(a) = — m / — cosh# ln(l + X(6\ a)) + 2V^ mcosh(yfc) 
L I 27t z ' 

J k= l 


The positions of the zeroes are fixed by the consistency conditions 


(24) 

(25) 


Y(xj±i— ;a) = — e la or — e 

X(y k±i^;a) = -1 


which follow immediately from the fact that (^3[) must still hold at the locations of 
the zeroes. We can write these relations as integral equations by using (G1.G3): 


K 


± a — (2 Mj + l)7r = i ln 
(26) 


fe=l 


tanh 

dO 


Xj - y k , 


tanh 


Xj + y fc mt 


1 


J 2ir sinh(xj — 0) 
~{2N k + 1)tt = mf?sinh(y fe ) 

J 


ln(l+X(0)) 


(27) 


+* 55 ln 

.7=1 
f dO 


tanh 

1 


yfc - Xj 


nr , 

— tanh 


y fc + Xj 77T 

2 + 4 


27T sinh(yfe — 0) 


In [(1 + e ia Y (0))(1 + e~ ia Y(0))] , 


where the Mj and N k are integers associated with a given solution. These integral 
equations are convenient for numerical solution. The branch of the logarithm is 
defined to be between — in and in; redefining this merely corresponds to redefining 
Mj and N k . Note that there are no poles in the integrands at 0 = Xj or 0 = y k 
because X(xj) = Y(y k ) = 0. We will show in the next section that continuity in 
a uniquely specifies the Mj and N k for a given a. Thus by studying how zeroes 
appear, one finds E(a) and hence 1(a) for all a. 

There is one special case of a where these equations are solvable in closed form. 
This is a = n, where T is Witten’s index. As discussed in |l^], we have X(0) = 0 
and Y(0) = 1 for all 0. This yields the required E\ = 0 as discussed in section 2. 
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4. 1(a) IN THE MASSLESS LIMIT 

In this section we will analyze the solutions of the TBA equations in the limit 
to —s- 0. We define 

A(0;a) = lim X(0 — ln(mI?/2); a) y(0\a) = lim Y(0 — 1ii(toI?/2); a) 

m — >0 m — >0 

As discussed in ||, these functions X and y are entire functions of 0 and a. In 
this section we exploit this analyticity to find integral equations which determine 
them for all a. This requires showing at which values of a zeroes enter the strip 
\ImO\ < 7t/ 2 and the corresponding Mj and N^. 

The functions X and y obey 


In A (6») 


(28) 

(29) In y(0) 


E i - / 0 — Xj 

intanh 1 


+ 


o=i 

r dff 


27t cosh(0 — O') 


In [(1 + e ia y(0'))(l + e~ ia y(0'))\ 


K 


Intanh 


fe=l 


Q - Vk 


dff 


27t cosh(0 — O') 


ln(l + X(0')) 


and 

K 

± a — (2 Mj + l)n = i Intanh 

k =1 

(30) + J Ssinh (lj-0) ln (l + ^)) 

—(2Nk + l)n = e Vk + i In tanh f ——— + — 

(31) + / ln [(1+e ” W))(1+e "“ y( ' ,))1 ■ 

Roughly speaking, these are “half” the solutions of the massive case. As mR —> 0, 
the functions X(0) and Y(0) are constant in the region |0| <C ln(2/ml?); this 
constant region becomes the asymptotic region 0 —> —oo in X(0) or y(0). 

Continuity requires that as we vary a, there are only two ways for a new zero 
to enter the strip \Im0\ < 7r/2. The first is for a single zero to enter at 0 = —oo. 
The second is for a pair of zeroes to enter, one at 0 = z + iw/2, and the other at 
0 = z — i7r/2. We adopt the conjecture of § that all zeroes in the strip occur at 
real values of exp(2(9/3). and therefore are on the real-0 axis. New ones must enter 
the strip at 0 = — oo. (The only exception to this statement is at a = 7r, where 
X{0) = 0 for all 0.) 

The first step to determining which zeroes are present at a given value of a is 
to solve the equations ( p8|]2g| ) in the limits 0 —> ±oo. This enables us to find the 
possible values of a where a zero enters or leaves at 0 = — oo. In subsequent steps 
we will show which of these places they do enter and leave (it turns out any time one 
can enter at —oo, it does). Since X(0) and y(0) approach constants as 0 —v — oo, 
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we can do the integrals in (|||) and (|fj|), yielding 

(32) [T(-oo)] 2 = 1 + 2cosa3 ; (-oo) + [3^(—cx))] 2 

(33) [^(-oo)] 2 = 1 + X(-oo) 

Requiring that the solution be positive at a = 0 hxes 

(34) X{-oo) = y(-oo) = 2cos(a/3). 

sm(a/o) 

In general, this means at a = (3 r + l)7r and a = (3 r + 2)i r, the number of zeroes 
(including those at oo) in X must change by l(mod2). Similarly, at a = (3r+3/2)7r, 
the number of zeroes in y must change by l(moc?2). For example, as a —> 97t/ 2, 
either an existing zero in y(9) returns to 9 —> — oo and goes away, or a new one 
enters. (We rule out pathologies like triple zeroes in the appendix.) 

Other useful pieces of information are the values as 9 —> oo, which are 

(35) X{9 —► oo) —> 2cos(a/2)e" exp(e) )V(oo) = 1 . 


The equation (g8|) only requires that that the coefficient in (|3|) be ±a/ 4 cos 2 («/2), 
but analyticity in a fixes the above result. For example, we see that a zero must 
appear in y{9) at a = 3 tt/2, since y(— oo) < 0 and y(9 —* oo) —> 1 for 3n/2 < 
a < 97t/ 2. Similarly, there must be a zero in X(9) when 2tt < a < 37t. Knowing 
the sign of the exponentially small term in (35) lets us see when zeroes go to +oo. 
Because the asymptotic behavior of X(9) at large 9 changes from 0^ to 0 T at 
a = 7T, 37 t , 57 t , ..., continuity requires that one of the zeroes Xj goes to oo and stays 
there. Once it is at oo, it no longer needs to obey (^). so it cannot move back to 
finite values. 

To find out what happens at higher values of a we must study the equations for 
the zeroes ( j30j ) and m in more detail. These equations simplify at the values of 
a where a zero is entering. For example, a zero yk enters at —oo, and in this limit 
the log terms become constants. Moreover, y{— oo) = 0 when a zero enters, so the 


integrand in ( |3l| ) vanishes as Xj —> —oo at these values of a. Thus when a zero yk 
enters or leaves at a = (3r + 3/2)7r by going to —oo, ( 0 ) simplifies to 

(36) (2 N k + 1) = J 


Because of the e Vk term in (30), we see that zeroes in y cannot go to +oo like those 
in X. 

Similarly, for a zero Xj entering (or potentially leaving) X at a = (3 r + l)7r or 
(3r + 2)7 t by going to —oo, (EhT) simplifies to 


(37) 


a — (2 Mj + l)7r = —7 tA' 


where the + sign in (^) must be chosen because as Xj increases from —oo, the 
terms on the right-hand side increase. (This can be checked by plugging in the 
expansion y(9) ~ 2 cos(a/3) + Aexp(20/3) + B exp(40/3) + ... valid for 9 —> —oo 
into the relation (f24|); A is given in (^) below.) We show in the appendix that 
choosing the + sign gives the correct contribution to the energy. For a zero Xj 
going to +00 at a = (2r + l)7r, one has 


(38) 


±a - ( 2 Mj + 1 )tt = 0 . 


Thus we see that when a zero Xj = — oo enters at a = a.j 1 the value Xj will go to 
oo at aj + Kjtt, where Kj is the number of zeroes in y at ctj. This phenomenon is 
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Location of first zero 



a/71 

Figure 1. The location of the first zero in X as a function of a 


displayed in fig. 1, where we plot the locations of various zeroes of X as functions 
of a. 

Now we can map out how at which values of a the zeroes enter at — oo, and at 
which values the zeroes in X approach +oo. For a < 7r, there are no zeroes, and the 
original TBA equations (H© apply. The situation at a = tt is atypical because 
of the special property that X(9) = 0 for all 9. Because K = 0 here, (|37]) and 
( j38| ) are the same, the zero with Mi = 0 enters at — oo and goes straight to +oo. 
The effect of a single zero x\ at infinity is to ensure that X(0) <0 for all 9. Thus 
the TBA equations for 7r < a < 37r/2 are almost the same as the ordinary TBA 
equations ( |I8| ) and (^9|): the only difference is that ln[A] is replaced by ln[— X], 

A zero y\ enters y at a = 37t/ 2. Because of the zero X\ already at infinity, J = 1 
here, and we have Ni =0. At a = 27t, X changes sign at —oo but not at +oo. 
Therefore, a zero Xi enters at — oo, and since K = 1 here, M 2 = 1. This zero then 
goes to 00 at a = 37 t, changing the ln[— X] back to ln[A]; effectively there is no 
zero in X for 3 t < a < 47 t. At a = 4w, a new one must appear because X(9) 
changes sign at 9 —> —00 but not near 00 . Still K = 1, so this new zero X 3 must 
have M 3 = 2. As we said above, at a = 97r/2, either a new zero j /2 appears, or the 
old one y\ goes away. The former must happen because J has changed from when 
yi entered; a zero with N\ = 0 can no longer satisfy (|3(]) and hence can not go 
back to — 00 . The new zero 1/2 has _/V 2 = 1. The zero X 3 with M 3 = 2 disappears at 
a = 5 tt, but again a new one must appear. This new zero has M 4 = 3, and since 
K = 2 here, it will go to infinity at a = 77 r. 

This pattern repeats modQTt. New zeroes enter X at a = (6r+l)7r, (6r+2)7r, (6r+ 
4)7 t and (6r + 5)7r with Mj = 4r, 4r + l,4r + 2 and 4r + 3, respectively. Zeroes enter 
y at a = (6r + 3/2) 7r and a = (6r + 9/2) 7r with Nk = 2r and 2r + 1, respectively. 
Zeroes Xj go to 00 at a = ( 6 r + l)7r, (6r + 3)7 t and (6r + 5 ) 71 . Thus when a goes 
up by 67 t, J increases by 4, K increases by 2, and the number of zeroes at infinity 
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e 


Figure 2. X(9) for several values of a 


goes up by 3. To summarize, 


while 


X2r—1 

enters at 

a = (3 r — 

2)7T 

X2r 

enters at 

a = (3j— 

1)7T 

Xj 

goes to oo at 

P 

II 

IS 

to. 

1)7T 

Vk 

enters at 

a = (3 k — 

3/2) 


Mj = j -1 
N k = k — 1 


All of this analysis is completely consistent with numerical solutions of the equations 
( p8p]p?M ) . We present some plots of the functions X(9) and y{9) for various 
values of a in figs. 2 and 3. 

Another very strong check of the results comes from comparison to perturbative 
results for 9 near — oo. In pi] . |2f| , it is shown that 

i [r(-|)] 2 r(f) 2e/3 

2 4/3 r(i + ^)r(i-£) 


(39) 


y{9) = 2cos(a/3) 
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e 


Figure 3. y(9) for several values of a 


(this corresponds to the coupling g = 2/3 in Q). This calculation is done by 
conformal perturbation theory and is completely independent of the TBA. In fact, it 
is not known how to extract this analytically from the TBA equations. However, we 
have checked that (^9|) is indeed satisfied by numerically solving the TBA equations 
for numerous values of a. 

This scenario is not the unique way of adding zeroes to the functions consistent 
with the relations ([57]|jf|[56|) , but it is the simplest. All the numerical evidence 
favors this scenario, as we have seen in the plots. In the appendix we calculate 
the energy, and show that for it to agree with the required value (10), many other 
possibilities can be ruled out. However, there seems to be a strange possibility that 
a zero could enter at +oo with a minus sign in ( psf ) . The analysis of the appendix 
does not rule out the possibility that at a = (6r + l)7r or (6r + 5)7r, additional zeroes 
simultaneously enter at +oo and — oo. While we cannot rigorously rule this out, 
we adopt the assumption of || that zeroes enter only from — oo, and our numerical 
work supports this assumption. 

These results are also completely in agreement with another conjecture of [|| 
concerning where zeroes enter their quantum monodromy operator. In our case, 
this agrees perfectly with our result that zeroes should enter y at a = (3r + 3/2)7 t. 











14 


PAUL FENDLEY 


5. Excited-state energies for arbitrary mass 

Once a mass is included, the functions X(9 ; a) and Y(9; a) remain analytic in 9 
for a given value of a. They also are continuous in to for fixed a. However, thay 
are no longer analytic in a, because they have discontinuities at certain values. We 
will utilize the results of the previous section to find which zeroes are present in 
the TBA equations with a mass ( |2~i| . p2| ) . In an appendix, we perform a similar 
computation for a free Dirac fermion. 

The functions X{9) and Y(ff) are even in 9. Thus zeroes appear in pairs (x^, — Xj) 
and (yfc, — y *.). If Xj satisfies (|2(]) with Mj . then —Xj solves it with —Mj. Because 
Yxm m ^ 0 X (6 — ln(2 /to); a) = X(9;a) and likewise for Y, the equations (|26|,|2^) 
smoothly turn into @01 as one takes the limit. Since the Mj and Nk are integers 
and the other pieces of (|2fi|. pTj) are continuous in Xj and ythe effect of turning 
on a very small mass is to shift the location of the zeroes and not change the Mj 
or Nk- So for m small enough, the zero structure is identical to that discussed in 
the previous section. Instead of a zero entering at — oo, we have a pair which enters 
at ±*7 t/2. As before, we discount the possibility that a quartet of zeroes enters at 
conjugate values (x + Itt/2,x — in/2,—x + in/2,—x — in/2) as mR is increased. We 
see no numerical evidence for a quartet entering. 

Our strategy in this section will be to fix a and start at mR small where we 
know the zeroes and the Mj and Nk- Once to is large enough, there is no guarantee 
that a given Xj or y k will remain real or even stay in the strip \Im(9)\ < n/2. For 
small mR , a zero of X(0) at 0 = x 3 corresponds to a zero of X{9) by the shift 
Xj ss Xj — ln(mf?./2). Thus the initial tendency as mR is increased is for the pair 
of zeroes to approach the origin. They can meet at the origin and move off along 
the imaginary axis. Eventually, the pair can reach 9 + iir/2 and 9 — iir/2 and thus 
not enter into the TBA equations any more. This is the only way they can go 
away: zeroes must appear in conjugate pairs, and we cannot have x r = x s because 
M r ^ M s when r ^ s (and likewise for y k )- Thus two zeroes cannot meet and go 
off into the complex plane except at the origin. 

If two zeroes meet at the origin, then equations ( pb| . p7| ) admit a solution Xj = 0 
or y j = 0. Because the integral vanishes when Xj = 0, and tanh(z+z7r/4) tanh(— z+ 
?'7t/ 4) = — 1, for Xj = 0 (|2tj) reduces to 

a — (2 Mfc + 1) 7r = —7 tK. 

This is identical to the condition ([n]) for a zero Xj to enter at —oo. At values of 
a where ( |37| ) is true, the zero isn’t there yet (it’s still at —oo). Thus there is no 
value of a where Xj = 0, and it is impossible for a zero in Xj to leave at any value 
of mR. Of course, (xj. —Xj) still moves to (oo, — oo) at the same values of a as in 
the massless case, because the equations for massive and massless are equivalent 
when 9 is large. 

It is possible for a zero y k to leave at some value of mR. If y k to be zero, then 

{2N k + 1) = J. 

This is the same as (^). However, as opposed to (|7|), this can be satisfied at values 
of a other than the value a k where the zero y k entered. As long as J hasn’t changed 
from its value at a k , a pair of zeroes can meet at the origin. Zeroes y k enter at the 
values a k = (3 k — 3/2)n. The next time the value of J changes is at a = (3 k — l)7r. 
Thus it is possible for zero y k to leave when (3 k — 3/2)7r < a < (3 k — l)7r. Our 
numerical evidence suggests that in this range it does in fact leave at some value 
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0 

Figure 4. Y(9) at a = 27r for several values of mR 


of mR. To illustrate this, in fig. 4 we present a plot of Y ( 9 ) for a = 27r for various 
values of mR. At mR ss .2, the zeroes meet at the origin and thus Y{6) >0 for all 
real 9. Then, for mR ~ .35, they leave the strip \Im6\ < 7r/2 and disappear from 
the TBA equations altogether. We see that Y(9) is continuous in mR for this fixed 
value of a. 

This observation lets us see at which values of a the functions are discontinuous 
when mR >0. At a = 27t, when mR > .4, there is no zero in Y(9). However, for 
a > 27t, there is a zero for any value of mR. This zero enters discontinuously, as 
we can see by looking at fig. 5, where y(9) is plotted at mR = 1 for a = 2n and 
2.00000017T. The functions are discontinuous even at mR < .2. There are similar 
discontinuities at the values a = (3 K — l)7r. 

Since the entire analysis of section 4 relied on continuity in a , let us reiterate how 
we found these discontinuities. In section 4, we utilized the fact that the functions 
are continuous in a in a special limit mR —► 0. This enabled us to find the which 
zeroes were present as a function of a. Then we fix a value of a. For a given a, the 
zeroes move around continuously as a function of mR, so the zero structure is the 
same for small enough values of mR, and we can thus exploit the mR —> 0 results. 
None of this precludes the possibility that if we fix mR and vary a, the functions 
can change discontinuously. In appendix B, we show that these discontinuities arise 
even for a free Dirac fermion. 

By examining the mR —> oo limit, we see that at a = (3K—1)tt, the energy shifts 
from that of a 2 (K — l)-particle state to a 2A"-particle state. The TBA equations 
are easily solvable in this limit, where X —► 0 and Y —> 1. The only contribution 
to the energy (|23|) comes directly from the zeroes yThe integral in (|3l]) goes to 
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0 


Figure 5. Y(9) at mR = 1 for a = 2ir and 2.00000017T 


zero in the large-mi? limit, and the mi?sinh(yfc) term must remain finite for there 
to be a solution. Thus yu —> 0 as 


(40) 


lim Yk 

mK—> oo 


(J — 2k + 1)7T 
mR 


where we used the fact that Nk = k — 1. In this limit the locations of the zeroes 
have the nice interpretation as the rapidities of the particles in the state. Because 
there are K zeroes for large mi? when (3A' — l) 7 r < a < (3 K + 2 ) 7 r, we have 


(41) lim E —> 2 Km. 

mR—* oo 

This corresponds to inserting a 2A'-particle state at the each end of the cylinder. As 
the cylinder radius R gets large, the particles are far apart so the energy approaches 
the constant value 2Km. The first correction is proportional to 1 /(mi? 2 ), as it 
should be for a collection of massive particles far apart: the momentum p ex 1/A, 
and for each particle E = •Jp 1 + m 2 ~ m + 0(l/(mA 2 )). All other terms in this 
expansion depend on the interactions between the particles; the full E gives the 
interaction energy non-perturbatively. One can compute the corrections perturba- 
tively by using methods described in j23|. In fig. 6 , we plot £3 as a function of mR. 


Let us examine £3 and £5 more closely. Both are two-particle states. There are 
two kinds of particle, the soliton and antisoliton. When rewritten in the appropriate 
N =2 language, soliton number is proportional to fermion number, and because they 
are a doublet under the supersymmetry, they must saturate the Bogomolny bound 
§. Since the states | n) have zero fermion number, the states |3) and |5) each consist 
of a soliton and an antisolition. From © we see that the individual momenta of 
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mR 

Figure 6 . E 3 at a = 37r as a function of mR. As mR gets large, 
E 3 « 2 yj 771 2 + ( n/R) 2 , the energy of two free particles of momenta 
±ir/R. 


the particles in the large R limit are ±n/R, in |3), and ±2n/R in |5). As discussed 
in j(|, one can presumably obtain the energy of two-particle states with higher 
momenta by using different values of Mj and N^. We have checked numerically 
that there are such solutions of the TBA equations, but have not pursued this line 
of investigation any further. 


6. Conclusion 

One application of this work is to the system of two-dimensional circular poly¬ 
mers. As discussed in [Q, partition functions calculated in this model are related to 
scaling functions for the number of closed self-avoiding random walks on a cylinder. 
The world-lines of the particles of this model correspond to the random walk, or 
equivalently the polymer. A single closed polymer will correspond to a two-particle 
state. Thus when we insert a two-particle state on each end of the cylinder, it 
corresponds to a circular polymer stretched from one end to the other. 

In this paper we have calculated various excited-state energies in the simplest 
integrable model with N =2 supersymmetry. We view this as a useful exercise, 
because it is a unitary theory, and because it shows that generalizing the work of 
II to models with non-diagonal S matrices is not prohibitively complicated. It 
seems likely that one could calculate the excited-state energies for any value of the 
coupling in sine-Gordon in this manner. Moreover, by choosing different values of 
the integers Mj and Nk, it seems likely that many more excited-state energies could 
be obtained. 

Our results also yield an interesting set of expectation values in N =2 supersym¬ 
metric theories which are independent of D-term perturbations like the “index” 
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studied in jlj] . That index had some remarkable properties in that the TBA equa¬ 
tions turned out to be equivalent to a differential equation, Painleve III. Moreover, 
because it depends on the parameters of the theory, it contains much information 
about the solitons of the theory. For example, it found useful application in four¬ 
dimensional gauge theories |2{|. This makes us hopeful that the results of this 
paper will also find future application. 


I would like to thank P. Dorey, K. Intriligator and H. Ooguri for helpful conver¬ 
sations. This paper was supported by NSF grant no. 9413208. 


Appendix A. The energy in the massless limit 


In this appendix, we calculate the energy in the m —► 0 limit and verify that it 
gives the correct answer © with our values of Mj and N k ■ This has the bonus of 
ruling out some alternative scenarios for the appearance of zeroes. 

We define £(a) = i?lim m _>o E(a). Then the expression (23) for the energy yields 


(42) 


£(<*) = — / e e L x (6) + 2j2< 


fc=l 


where L x {9) = ln(l + X(9)). Similarly, we define Ly{9) = ln[(l + e la y(9))( 1 + 
e~ la y{9))\. To simplify this, we use a well-known TBA trick. We take the derivative 
with respect to 9 of (|28|) and use this to eliminate the e 6 in ([42|) , yielding 


K 


£ = 2 E' 


+ 


k =1 


d9 5X(9) 1 

~ 89 X{9) 


j 


£*(*)-E 

i=i' 

d9 d9' 


d9 


1 


-L x (9) 


7r sinh(0 — Xj ) 

1 SL X (9') 


7 r 27 t X ^ ^cosh(0 — 9') 


69' 


The double integral in the final term can be simplified back to a single integral by 
using (E9h. After integrating by parts, one finds 


£ = 


Vk 


d9 6X{9) 1 


-L x {9) 


d9 


K 


K 

2 E e 

k =1 

J 

-E 

3 =1 

-iln[(A(-oo)) 2 ]ln[((-l)^(-oo)] 


7r sinh(0 — Xj) 


69 X(0) 

r-7T L 'V(e) - £ 


fc =1 • 


f d9 6y{9) 1 
7T 89 y{9) 

d9 1 

7T sinh(0 - y k ) 


Ly(9) 

Ly{9) 


where we simplified the final term (the surface term from the integration by parts) 
by using (^2 33). The zeroes at infinity (there are Joo of them) do not contribute, 
so we can rewrite the sum over j to start at + 1. Notice we can change variables 
to rewrite the first two integrals as simple integrals, so for example 

f d9 5X(9) 1 r f x{oo) du 1 , . 

/- " y(o\ Lx ^ = / -ln(l + u). 

J TT 69 X(9) J X (- 00 ) 7T u 
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We use the equations for zeroes (|3(|j31| ) to eliminate the terms with the 1 / sinh and 
the e Vk , yielding 

r x (°°) rjn, i ry{°°) d v 1 

£ = / -ln(l +u) + / — -ln\{l + e ia v)(l+ e - ia v)} 

Jx(- oo) u Jy(-oo) 7T v 


-- ln[(*(-oo)) 2 ] ln[y(-oo)(-l) K ] 

7T 

K 


) K ) 

j 


(43) 


+2JKtt — 2 + l)7r + 2 ^ (±a - {2Mj + 1)tt) 

fe=1 j = Joo+l 

Because the first three terms are always finite and depend only on 3f(±oo) and 
(V(±oo), they are periodic under in a —> a + 67 t. However, as discussed in section 
2 , we know that 

(44) 

and the full answer is not periodic. This means that zeroes must enter into X and 
y in order so that the last three terms make £ (a) continuous in a. This puts strong 
constraints on the values of Mj and TV* and at which values of a the zeroes enter. 
Despite the fact that it looks like (]43|) will be discontinuous when a zero enters, 
the conditions (^,|^,| 8 |) ensure continuity. This condition of continuity rules out 
many other conceivable ways zeroes could enter. 

To simplify it further, we take the derivative of £ (a) with respect to a. By using 
the relations (|||) and (^) and their derivatives with respect to a, most of the 
resulting terms cancel, leaving only 


6 £ 

5a 


dv 


2 sin a 


/2 cos(a/3) 7T 1 + 2 COS av + V 2 


J 

E 

j = J a 0+1 


(± 1 ). 


For a < 7r, one can do the integral, yielding —a/Sir in agreement with ([44|). However 
at integer values of a/ 7 r, one must carefully check if there are discontinuities because 
the denominator and numerator both can become zero. One sees that at a = 
( 6 r + 2 ) 7 t and ( 6 r + 4)n, the integrand is 2 5(v + 1) giving discontinuities of —1 
in 5£/5a. Similarly, at a = ( 6 r + 3 ) 7 t there are discontinuities of +1. These 
discontinuities must be cancelled by the second term. This means a zero with +a 
in (^o|) must enter at a = ( 6 r + 2 ) 7 r and ( 6 r + 4 ) 7 T, while at a = ( 6 ?' + 3 ) 7 r, either one 
with —a should enter, or one with +a should go to oo, where it does not contribute 
to (|43[). Since we show in section 4 that a single zero cannot enter at a = ( 6 r + 3 ) 7 r, 
the latter possibility must hold. 

These conditions are satisfied by the scenario discussed in section 4, and we take 
it as strong evidence that it is correct. The analysis of this appendix rules out, for 
example, the possibility that a triple zero enters at a = (6 r + l)7r. 


Appendix B. The excited states of a free Dirac fermion 

In this appendix we derive the excited-state energies of a free Dirac fermion by 
applying the general methods of this paper. This will clarify several issues which 
arose in the interacting case, in particular the nature of the discontinuity discussed 
in section 5, where the energy is discontinuous at particular values of a when the 


mass is non-zero. 
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Just like the interacting case, the expectation value of interest is (||). This cor¬ 
responds to placing the fermion and antifermion at opposite (imaginary) chemical 
potentials. Although (f]) looks periodic in a, we define the values for a > n as 
the appropriate continuation of those for a < n. Since the particles are free, we 
can evaluate the trace explicitly. We define a function T(9) which satisfies the 
functional relation 

(45) T(9 + %n/2)T{9 - in/2) = e - m - Rcoshe + 2 cos a + e mRcoshf> . 


It follows that T{9 + 2n) = T(Q). In this non-interacting case, we can find all the 
zeroes of T explicitly. The relation (fit]) requires that 

exp (ia + mR cosh(t„ + in/2)) = —1, 

where T(t±) = 0 and T(—t±) = 0. Taking the logarithm, we have 

(46) a + mR sinh t+ = —(2n — l)n 

(47) — a + mR sinh t~ = — (2n — l)n 


where n is some integer. Using the Lemma of section 3, one can write an integral 
equation for InT: 


ln[T( 0 )] 

(48) 


mR cosh 9 + In 


tanh 


9 - t fc 


tanh 


9 + tfc 


d9' 


1 


2n cosh(0 — 9') 


In [(i + e i« e -™*™she 


_|_ g—ioig —mR cosh 6^ 


where the sum over k is only over those zeroes which are in the strip —n /2 < 
Imtk < n/ 2 . 

The energy levels are found by the large-0 asymptotic expansion as in © 1 - A 
bulk term is present here which does not appear in the supersymmetric case. This 
term is independent of a, so we won’t bother to discuss this subtlety. The energy 
is 


(49) E(a) = Eo(a) +2 m cosh(tfc) 

k 

E 0 {a) = - m /fcosh01n [( l + e^—)(l + e^e- co -)] 

where again the sum over k is over only the zeroes in the strip. Notice that Eq{q) 
is periodic in a, and denotes the usual Casimir energy of a Dirac fermion. 

We now need to decide which zeroes are present for a given value of a. At a = 0, 
Eo is indeed the correct ground state energy for a Dirac fermion, so there are no 
zeroes in the strip here. Thus the zeroes t^ must solve ( ptilfl’ij ) with n > 1, so 
that Imt^ = n when a = 0. (The relation ( fi5| ) only requres that T{ t^) = 0 or 
T(t^ + in) = 0, so it is consistent to require that T(t^) = 0 only for n positive. 
One can find different excited states of the system by restricting the allowed n to 
be a different subset of the integers.) 

As discussed in ||, in the massless limit T(9;a) = lirn m ^ 0 ^(^ ~ ln(mi?/2); a) 
is analytic in 9 and a. The zeroes of T must therefore obey 

a + e* n = —(2 n + l)7r 
—a + e* n = —{2n + l)7r 
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with n > 1 for all a. For a < ir, all the have imaginary part in, and do not 
contribute to the energy (|50|). At a = n, = — oo. As we go past a = n, this 
zero “jumps” onto the real axis. Even though this seems discontinuous, one can 
easily check that T is continuous: the extra term coming from the zero in (ff 8 | ) is 
zero when tk = — oo. Similarly, at a = (2 K — l) 7 r, the zero jumps onto the real 
axis and hence enters into expressions for T and the energy. If we were to take a 
to —( 2 K — l) 7 r instead, then the zero tjr jumps onto the real axis. 

A similar phenomena happens in the massive case. At a = (2 K — l)n, we see 
from ([id]) that the zero t j". jumps from i 7 r to 0. This solution reduces to the massless 
one appropriately in the limit m —► 0. (It is possible to find a continuous solution 
to (]f7|) by having the zeroes and —tj, just switch places at a = (2 1\ — l)7r, but 
this solution is periodic in a and does not reduce to the correct massless solution.) 
Thus for a > (2 K — l)7r, the zero tj(- must now must be included in the sum in (48) 
and likewise in ( |50| ) . However, the interesting distinction between the massless and 
massive cases is that T and E are not continuous at these values of a, as is easily 
seen by examination of ( fl 8 | ) and (50). 

The physical interpretation of this is clear. For {2K — l)n < a < (2 K + l)-7r, the 
energy is that of a state of K fermions and K antifermions. The zeroes ±t& with 
k < K are the rapidities of the particles, so that the energy of a given particle is 
mcoshtfc. The energy still contains the Casimir piece Eo(a); because the particles 
are non-interacting it is not changed by presence of physical particles in the state. 
The example treated in the main part of this paper is interacting, so the energy 
only decouples into a “Casimir” piece and an “excited-state” piece in the large-1? 
(dilute) limit. However, we have seen in this appendix that the discontinuities in 
the energy and in the function used to determine it are very natural physically, if 
somewhat surprising mathematically. 
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